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Abstract. Gear box dynamics is characterised by a periodically changing stiffness. 
In real gear systems, a backlash also exists that can lead to a loss in contact be- 
tween the teeth. Due to this loss of contact the gear has piecewise linear stiffness 
characteristics, and the gears can vibrate regularly and chaotically. In this paper we 
examine the effect of tooth shape imperfections and defects. Using standard methods 
^""^ ' for nonlinear systems we examine the dynamics of gear systems with various faults 

, in meshing stiffness. 
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Gears are very common systems, and practically impossible to replace 
in various applications where mechanical power must be transferred. 
Time varying mesh stiffness due to multiple teeth contact and a back- 
lash between the teeth give rise to complex behaviour [1-6]. In conse- 
quence, under a dynamic load, a typical gear system is a nonlinear 
oscillator, exhibiting a range of complex behaviour including chaos 
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^ . [4,7-13]. In use the geometric parameters of the gears change, and 

this causes the corresponding nonlinear response to change [4-5,14-15]. 
Choy et al. 16 and Kuang and Lin |17j examined the effect of tooth 
wear. Vibrations have also been modelled by including stochastic forces 
[4,14-15,18]. 

In practice it is important to minimise the effect of noise and keep 
the machine as close as possible to a stable response. In this paper we 
classify meshing faults and examine the effect of broken teeth and mesh- 
ing stiffness fluctuations on the vibration response. The possiblility of 
amplitude jumps in systems with meshing defects is demonstrated. 

© 2008 Kluwer Academic Publishers. Printed in the Netherlands. 
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Figure 1. One stage gear system. 

2. Modelling of gear dynamics 

Consider the single gear-pair system shown in Fig. 1. In non-dimensional 
form, the equation of motion can be written [3J Ej as 



d 2 2C d k(r)g(x, V ) 
x H — x H t, 



B + Bxcos{t + 6) 



dr 2 



u) dr 



where 
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riMi/J 1 -r-r 2 M2/J2. 



uj, C) ^( r )) 9( x iV)i V an d -B(t) ^ and the other symbols are defined in 
Table 1 and shown on Fig. 1. 

In the analysis that follows, the stiffness functions k(r) and g(x,rj) 
need special attention. g(x, rj) has a piece-wise character due to the 
backlash rj, and is shown in Fig. 2. k(r) is the meshing stiffness arising 
from the interaction of a single-pair or multiple teeth in contact. For an 
ideal gear system we have followed references (U EI] and assumed that 
this meshing stiffness changes periodically. Possible variations from the 
ideal case, and other possible meshing errors, are plotted in Fig. 3. 

Figure 4 shows the results of simulations of the model given by 
Eq. (1), with time dependent meshing stiffness but without errors. We 



lit_fri_7.tex; 8/02/2008; 19:54; p. 2 



-20 -15 -10 -5 5 10 

x 

Figure 2. The nonlinear stiffness function g(x,ri). 

Table I. Symbols and parameters used in the analysis. 



h, h 


moments of inertia 


Ipl, 1p2 


rotational angles 


x — ritpi — r-2ip2 


relative displacement 


V 


relative velocity 


Xo, v 


initial conditions 


Mi, Mi 


external torques 


UJ 


excitation frequency 


T 


dimensionless time 


c 


damping 


V 


backlash 


k(r) 


meshing stiffness 




nonlinear stiffness function 


B, Bo, B\ 


external excitation 




distance between incerasing teeth contacts 




standard deviations 



have used following system parameters: u = 1.5, ( = 0.08, Bo = 1.0, 
B\ = 4.0, r] = 10. With any nonlinear systems multiple solutions may 
coexist, and the solution obtained depends on the initial conditions. 
With the above parameters, there are indeed multiple solutions for the 
gear model, and this effect was examined in detail in a previous paper 
|19j . In Figs. 4a and b we show regular and chaotic solutions, depending 
on the initial conditions for [xo,vo] = [—9,1] and [xo,fo] = [— 9, — 1], 
respectively. Figure 4c shows the two coexisting attractors, obtained 
for various initial conditions, on the same graph. 
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Figure 3. Various realizations of meshing stiffness k(r). '1' corresponds to the ideal 
system without errors while '2' and '3' show the meshing stiffness with one and two 
broken teeth. '4' has a randomised distance Si between increasing teeth contacts with 
a standard deviation as = 0.2<5; (Si — 0.8n in the non-dimensional time domain), and 
'5' is a randomised meshing stiffness. Here the amplitude changes with a standard 
deviation a\ = O.f related to the maximum deterministic value k max = 1. 



3. Errors in meshing stiffness 



In this section we examine the effect of meshing stiffness errors. First 
consider a gear with one or two next neighbour teeth missing on one 
of the gear wheels, where each gear wheel has 50 teeth. The shape 
meshing stiffness fc(r) are given in Fig. 3 as '2' and '3', respectively, 
comparing to the ideal case '1'. The response is simulated using the 
model given in Eq. 1 and modelling the meshing stiffness fc(r) by curves 
'2' and '3' (Fig. 3). Although the effect of one broken tooth appears to 
be fairly benign in terms of the gears dynamics, if two next neigbour 
teeth are broken the result is a complex response of the system showing 
the characteristic amplitude jump phenomenon as the solution changes 
from the regular to the chaotic attr actor. This is visible in Fig. 5, which 
shows a time history for this case. Also shown is the reverse jump from 
the regular to the chaotic attractor, but it is clear that the system 
stays for longer time in the chaotic attractor with intermittent regular 
motion. This result confirms previous results on stochastic jumps |14l H| 
in systems with a stochastic force. However, the system modelled here 
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Figure 4- Time series and the Poincare map for an ideal system for various initial 
conditions [a;o,t>o]; regular motion for [xo,«o] = [—9, 1] (a) and chaotic motion for 
[xo,t>o] = [ — 9j ~1] (b), respectively. Poincare map for various initial conditions (c). 
The characters a and b in figure c denote the regular and chaotic attractors give by 
time series (a) and (b). 



is fully deterministic and the broken teeth act as additional parametric 
excitation. 

We have also investigated vibrations of gears with a random distance 
5i between their increasing teeth contacts (Fig. 3, '4'). The results for 
two different noise levels and two different initial conditions [xo,^o], 
which correspond to different attractors in deterministic case (Fig. 4), 
are shown in Figs. 6a-d. Interestingly, for weak noise the system chooses 
the chaotic attractor. This conclusion differs from that obtained in the 
paper by Warmihski et al. jlj but the assumptions about the noise 
are different. Warmihski et al. [I] used an external stochastic force 
generated by stochastic Langevin simulations rather than the stochastic 
stiffness modelling in the present paper. For stronger noise the mo- 
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Figure 5. Time series for a gear system with two broken neighbouring teeth. Initial 
conditions: [a;o,«o] = [—9, 1]. 
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Figure 6. Time series for a gear system with a randomised distance <5; between 
increasing teeth contacts with two different standard deviations: as = 0.25 in (a) and 
(b), and as = 0.35 in (c) and (d); and two different initial conditions [xo, vo] = [—9, 1] 
in (a) and (c), and [a;o,«o] = [—9, —1] in (b) and (d). 
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Figure 7. Time series for a system with randomised amplitude showing regu- 
lar motion for initial conditions [a;o,Uo] = [—9,1], (a), and chaotic motion for 
[x ,vo] = [-9,-1], (b). 
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Figure 8. Bifurcation diagram for the case of randomised meshing stiffness am- 
plitude: x is plotted stroboscopically against the square deviation of fluctuating 
meshing stiffness au for two different initial conditions [xo,uo] = [—9, 1] (the upper 
panel) and [a?o,«o] = [— 9,— f] (the lower panel). 



tion shows an intermittent character with short jumps into the regular 
attractor as in the previous case with broken teeth (Fig. 5). 

Figure 7 corresponds to the meshing stiffness with randomized am- 
plitude but regular distance (Fig. 3, '5'), for = 0.1. Here we observe 
that neither attractor is favored for these conditions. Eventually, for 
a sufficiently large noise level (cr^ « 0.3) the response returns to the 
intermittent behaviour with jumps between the regular and chaotic 
attractors. This effect is shown clearly in Fig. 8, which shows the 
bifurcation diagram x modulo 2tt against noise level a^. 



4. Conclusions 



We have examined the dynamics of gears in the presence of meshing 
faults. Such faults may arise due to wear during or incorrect tolerance in 
their production. The analysis of various types of errors and tooth faults 
highlights the presence of dynamic jumping phenomenon. Such jumps 
between different the types of motions, namely chaotic and regular, 
can be crucial for the system reliability. In this respect our results are 
consistent with earlier results [1411181 H|. Moreover, the system is more 
sensitive to errors in the distance between teeth than fluctuations in 
the stiffness magnitude, although the qualitative effect is similar. One 
broken tooth has little influence on the dynamics of thegears, although 
two broken teeth can have a significant effect. 
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